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Measurement of the A° lifetime using semileptonic decays 
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We report a measurement of the Aj lifetime using a sample corresponding to 1.3 fb~^ of data 
collected by the DO experiment in 2002-2006 during Run II of the Fermilab Tevatron collider. The 
A° baryon is reconstructed via the decay A^ fii>A'^ X. Using 4437 ± 329 signal candidates, we 
measure the A° lifetime to be t(Aj) = 1.290 Iq iio (stat) lo ogi (syst) ps, which is among the most 
precise measurements in semileptonic A° decays. This result is in good agreement with the world 
average value. 

PACS numbers: 14.20.Mr, 14.40.Nd, 13.30.Eg, 13.25.Hw 



Lifetimes of b hadrons provide an important test of 
models describing quark interaction within bound states. 
The experimental measurement of the lifetimes are in 
reasonable agreement with the theoretical predictions fT, 
[H, 3], but further improvement in the experimental and 
theoretical precision is essential for the development of 
non-perturbative quantum chromodynamics. 

The lifetime of b baryons recently attracted a spe- 
cial interest. The current world average A^ lifetime 
is r(A°) = 1.230 ± 0.074 ps, and the ratio of the 
baryon and i?" meson lifetimes is r(A°)/T(i?°) — 0.80 ± 
0.05 3 1 ill good agreement with the theoretical predic- 
tion r(Ag)/T(B°) = 0.86 ± 0.05 [5]. However, the recent 
A[J lifetime measurement from the CDF collaboration in 
the A^ — > J/i/jA decay gives a significantly larger value: 
T{A"^,) = 1.593toml ± 0-033 ps [5], not included in the 
quoted world average. Additional AJJ lifetime measure- 
ments could provide a potential resolution of this incon- 
sistency. 

This Letter presents a measurement of the A° lifetime 
using the semileptonic decay A° ^PA+AT, where X is 
any other particle. Charge conjugated states are implied 
throughout this paper. The A+ baryon is selected in the 
decay A+ — > Kgp. The sample corresponds to approxi- 
mately 1.3 fb~^ of data collected by the DO experiment 
in Run II of the Fermilab Tevatron Collider. 

The DO detector is described in detail elsewhere Q. 
The components most important to this analysis are the 
central tracking and muon systems. The central track- 
ing system consists of a silicon microstrip tracker and 
a central fiber tracker, both located within a 2 T super- 
conducting solenoidal magnet, with designs optimized for 
tracking and vertexing at pseudorapidities \ri\ < 3 and 
\ri\ < 2.5 respectively (where 77 = — ln[tan(0/2)] and 9 is 
the polar angle of the particle with respect to the proton 
beam direction). The muon system is located outside the 
calorimeters and has pseudorapidity coverage \r]\ < 2. It 
consists of a layer of tracking detectors and scintillation 
trigger counters in front of 1.8 T iron toroids, followed by 
two similar layers after the toroids Q . The trigger system 
identifies events of interest in a high-luminosity environ- 
ment based on muon identification and charged tracking. 
Some triggers require a large impact parameter for the 
muon. Since this condition biases the lifetime measure- 
ment, the events selected exclusively by these triggers are 
removed from our sample. All processes and decays re- 
quired for this analysis are simulated using the evtgen 



[sl generator interfaced to PYTHIA [oj and followed by full 
modeling of the detector response using GEANT [lOj and 
event reconstruction. 

Reconstruction of the A° decay starts from the selec- 
tion of a muon, which must have at least two track seg- 
ments in the muon chambers associated with a central 
track, with transverse momentum px > 2.0 GeV/c. All 
charged particles in the event are clustered into jets using 
the Durham clustering algorithm [ll|. The products of 
the A+ decay are then searched for among tracks belong- 
ing to the jet containing the identified muon. 

The primary vertex is determined using the method 



described in Ref. I4I. The Kg meson is reconstructed 



as a combination of two oppositely charged tracks that 
have a common vertex displaced from the pp interaction 
point by at least four standard deviations of the measured 
decay length in the plane perpendicular to the beam di- 
rection. Both tracks are assigned the pion mass and the 
mass of the tt^'^tt^ system is required to be consistent with 
the Kg mass to within 1.8 standard deviations. Combi- 
nations consistent with the A — > pn hypothesis, when 
either track is assigned the proton mass and the mass 
of the pir system lies between 1.109 and 1.120 GeV/c^, 
are rejected. Any other charged track in the jet with 
Pt > 1.0 GeV/c and at least two hits in the silicon detec- 
tor is assigned the proton mass and combined with the 
neutral extrapolated Kg candidate to form a A+ can- 
didate. Their common vertex is required to have a fit 
X^/d.o.f.< 9/1. The A+ candidate is combined with the 
muon to make a Ajj candidate, and its invariant mass 
is required to be between 3.4 and 5.4 GeV/c^. A com- 
mon vertex for the A+ candidate and muon is required 
to have a fit x^/d.o.f. < 9/1. The transverse distance 
between the A° and A+ vertices is calculated and is 
assigned a positive sign if the A^ vertex is closer to the 
primary vertex, and a negative sign otherwise. The A° 
candidate is required to have —3.0 < d^/a{d^) < 3.3, 
where (j{d^r^) is the uncertainty of the d^ measurement. 
The upper bound on the distance between A° and A+ 
vertices reduces the background significantly, since the 
A^ lifetime is known to be very small: 0.200 ± 0.006 
ps 4]. 

To further improve the A° signal selection, a likelihood 
ratio method [13j is utilized. This method provides a sim- 
ple way to combine many discriminating variables into a 
single variable with an increased power to separate signal 
and background. The variables chosen for this analysis 
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are the A° isolation, the transverse momentum of the Kg, 
proton and Aji" candidates, and the mass of the /^A+ sys- 
tem. The isolation is defined as the fraction of the total 
momentum of charged particles within a cone around the 
/iA+ direction carried by the A° candidate. The cone is 
defined by the condition ^ (Atj)^ -f (A^)^ < 0.5, where 
Ar] and Acj) are the difference in pseudorapidity and az- 
imuthal angle from the direction of the A° candidate. 

Figure [1] shows the invariant mass M{Kgp) for the 
selected A° candidates. The fit to this distribution is 
performed with a signal Gaussian function and a fourth- 
order polynomial function for the background. The A^ 
signal contains 4437 ± 329 (stat) events at a central mass 
of 2285.8 ± 1.7 MeV/c^ The width of the mass peak is 
a — 20.6 ± 1.7 MeV/c^ consistent with that observed in 
the simulation. 

Simulation shows that the contribution from the i?^ — > 
KgTT decay when a pion is assigned the proton mass has 
a broad M{Kgp) distribution with no excess in the A+ 
mass region. 
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FIG. 1: The ifgp invariant mass for the selected Aj^ candi- 
dates and fit overlaid (see text). Notice the suppressed-zero 
scale of the vertical axis. 



Since the final state is not fully reconstructed, the A° 
proper decay length cannot be determined. Instead, a 
measured visible proper decay length A^'^ is computed as 
A*^ = mc [Lt ■ pT(/i^?)) / Ipt(M? )P- J^T is the vector 
from the primary vertex to the h.^ vertex in the plane 
perpendicular to the beams, pt(mA+) is the transverse 
momentum of the /xA+ system and m = 5.624 GeV/c^ is 
taken as the A° mass 0|. 

To determine the A" lifetime, the selected sample is 
split into a number of A*^ bins. The mass distribution 
in each bin is fitted with a signal Gaussian and a fourth 
degree polynomial background. The position and width 
of the Gaussian are fixed to the values obtained from 
the fit of the entire sample (see Fig. [1]). The Gaussian 
normalization and background parameters are allowed to 
float in the fit. The range of A*^ and the number of signal 



events fitted in each bin together with its statistical 
uncertainty Ui are shown in Table HI 



TABLE I: Fitted signal yield in different A^'' bins 
A'*^ range(cm) Number of signal candidates rii ± Oi (stat) 



-0.06, -0.04] 


62 ±48 


-0.04, -0.02] 


66 ±69 


[-0.02, 0.00] 


587 ±156 


[0.00, 0.02] 


1172 ± 173 


[0.02, 0.04] 


999 ± 99 


[0.04, 0.06] 


540 ± 69 


[0.06, 0.08] 


299 ± 54 


[0.08, 0.10] 


225 ± 44 


[0.10,0.20] 


454 ± 64 


[0.20, 0.30] 


47 ±34 



The expected number of signal events in each bin n1 
is given by n1 = Ntot /■ /(A*^)dA*^, where Ntot is the 
total number of /iA+ events, and f{\^) is the probability 
density function {pdf) for A^^. The integration is done 
within the range of a given bin. 

In addition to AJJ fiDA'^X decays, the A+ baryon 
can also be created in cc or 66 production, along with a 
muon from the decay of the second c or 6 hadron. In what 
follows, these processes are referred to as peaking back- 
ground, since they produce a A^ peak in the Kgp mass 
spectrum imitating the signal. Such events are recon- 
structed as Aj] candidates, and have a fake vertex formed 
by the intersection of the muon and A+ trajectories. The 
simulation shows that the distribution of A*^ for such a 
fake vertex has a mean of zero and a standard deviation 
of sal50 fiui. 

The expression for /(A*^) takes into account 
the contribution of signal and peaking background: 
/(A^^) = (1 - rbck)/sig(A^O + r-bck/bck(A*0- Here rtck 
is the fraction of peaking background, and /sig(A*^) and 
/bck(A*''^) are the pdfs for signal and background respec- 
tively. The background pdf is taken from the simulation. 
The signal pdf is expressed as the convolution of the de- 
cay probability and the detector resolution: /sig(A*^) = 
/ dKH{K) [diX)K/{cT) exp(-XA/(cT)) (g> R{\^' - A, s)] . 
Here, r is the A^ lifetime, and 6{X) is the step function. 
The factor K = pT{nh.c) I Pt{J^V} ^ measure of the 
difference between the measured pt(a'A+) and true 
momentum of the A° candidate, and H{K) is its pdf. 
The i?(A*^ — A, s) is a function modeling the detector 
resolution. A scale factor s accounts for the difference 
between the expected and actual A*^ resolution. 

The H{K) distribution is obtained from the simu- 
lation. The contribution of decays l^-^^t 

fiDT^cTT with Sc ~> A+TT is taken into account. The 



A+fi*' with the decaying 



contributions of A^ 
semileptonically, — > fiPAcX and Aj] t^DA'^ with 
T~ — > yb~i>^VT are found to be strongly suppressed by the 
branching fractions and low reconstruction efficiency. To 
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obtain H{K), the K factor distribution of each process is 
weighted with its expected fraction in the selected sam- 
ple. This is computed taking into account both the re- 
construction efficiency and the branching fraction of each 
process. The fraction of £~DA'^ in semileptonic de- 
cays has been measured recently to be 0.47lo;}o [I]- We 
use this result in our analysis. 

The resolution function is given by i?(A*^ — A, s) = 
/ /res(o')G(A^^ — A, CT, s)da, where /res(o') is the pdffoi the 
expected resolution of A*'^, and G is a Gaussian function 
G{X^^ - A, a, s) = l/iV2^as) exp[-(AAf - X)yi2ah^)]. 
The (Ts is the decay length uncertainty, which is deter- 
mined for each candidate from the track parameter un- 
certainties propagated to the vertex uncertainties. 

To determine /resCf), signal and background subsam- 
ples are defined according to the mass of the Kgp sys- 
tem. All events with 2244.7 < M{K^p) < 2326.9 
MeV/c^ are included in the signal subsample, and all 
events with 2183.9 < M{K'^p) < 2225.0 MeV/c^ and 
2346.6 < M{KIp) < 2387.7 MeV/c^ are included in 
the background subsample. In addition, the events in 
both subsamples are required to have a measured proper 
decay length exceeding 200 /im. This cut reduces the 
background under the A+ signal and the contribution of 
peaking background. The fics{<^) distribution is obtained 
by subtracting the distribution of expected resolution in 
the background subsample from the distribution in the 
signal subsample. 

The A^ lifetime is determined by the minimization of 
— J2i'''''"' where the sum is taken over all 

bins of measured proper decay length (Table [T|. The free 
parameters of the fit are Ntot, ■''(A°) and rbck- A separate 
study is performed to measure the resolution scale factor 
using the decay D*+ D°tt+ with D° fi+vK^i:-. It 
has a similar topology to that of the A° — > fiPA'^ decay. 
Since the D*~^ meson comes mainly from cc production, 
its decay vertex coincides with the primary interaction 
point. The distribution of the D*^ proper decay length 
is mainly determined by the detector resolution and can 
be used to measure the resolution scale factor. A value 
of 1.19 ± 0.06 is found. The scale factor in the lifetime 
fit is fixed to this value and varied later in a wide range 
to estimate an associated systematic uncertainty. 

The lifetime fit gives r(A°) = 1.290tEJ "i^ (stat) ps, and 
the fraction of peaking background rbck = 0.160^0 q74 
(stat). Figure [2] shows the distribution of the number 
of A+/i events versus A^^ together with the result of the 
lifetime fit superimposed. The lifetime model agrees well 
with data with a x^/d.o.f.= 5.5/7. The dashed line shows 
separately the contribution of the peaking background. 

The method used to fit the mass distribution in each 
of the A*'^ bins is the most significant source of system- 
atic uncertainty. The fit sensitivity is tested by refitting 
each A*^ bin for the mass interval between 2.17 and 2.40 
GeV/c^ with a linear parametrization of the background. 
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FIG. 2: Measured fiA^ yields in the A bins (points) and 
the result of the lifetime fit (solid histogram). The dashed 
histogram shows the contribution of peaking background. 



Binning effects of the mass histograms are checked by 
performing fits to the data with bins of half the nomi- 
nal width and with the lowest and highest bins excluded. 
The lifetime fit is performed again for each test. The 
largest deviation of r(A°) is 0.067 ps, which is given as 
the systematic uncertainty due to the mass-fitting proce- 
dure. The parameters describing the peaking background 
are varied by their uncertainties. The largest shift in the 
fitted lifetime is 0.012 ps. 

The selected sample can also contain a contribution 
from B — > iiDA'^X decay. Its branching fraction is 
unknown; only the upper limit Br(i3 — > eDA'^X) < 
3.2 X IQ-^ at 90% CL is available The possible con- 
tamination from this decay would reduce the fitted A^ 
lifetime, since the K factor for these events is smaller. 
The upper 90% CL limit on the fraction of this decay in 
the selected sample is estimated to be 5%, which would 
result in the reduction of the A^ lifetime by 0.027 ps. 

The value of the scale factor is varied by ±20%, and 
shifts of approximately ±0.036 ps are observed in the fit- 
ted lifetime. This value is also included in the systematic 
uncertainty. 

The fraction of A° M^^^ decay in the semileptonic 
A^ decays is varied between 0.3 and 0.6. The lower bound 
is selected to be larger than the current uncertainty in 
this fraction Q to take into account the possible contri- 
bution from decays to tvA'^ and other heavier states with 
lower mean K factor. The shift of 0.025 ps in the fitted 
lifetime is taken as the systematic uncertainty due to the 
branching fractions in the K factor. The mean of the K 
factor distribution does not change significantly with the 
Pt of the muon, however the shape of the distribution is 
changed. To estimate the possible variation of the A^ life- 
time, the distribution for ^iDA'^ decays is generated with 
a cut of pri^J-) > 6 GeV/c and the fit is repeated. A shift 
of 0.005 ps is observed, which is assumed as the uncer- 
tainty due to the momentum dependence of the K factor. 



The change in the K factor distribution due to the un- 
certainty in generation and decay of B hadrons has been 
estimated in other analyzes to be less than 2% 14, 1^. 
Therefore we shift all K factor values by ±2%, and ob- 
serve a shift of 0.026 ps in the fitted lifetime. The overall 
systematic uncertainty due to the K factor distribution 
is estimated to be 0.036 ps. The effect on lifetime mea- 
surement due to misalignment of elements of the tracking 
detector is determined by rescaling the geometrical posi- 
tion of all detectors within uncertainties of the alignment 
procedure. The resulting variation of the A|J lifetime is 
estimated to be 0.018 ps. 

The systematic uncertainties are summarized and 
added in quadrature in Table [ill Total systematic un- 
certainty of this measurement is estimated to be 0.09 ps. 

In addition, several consistency checks of this analy- 
sis are performed. The fitting procedure is applied to 
the simulated A° f^i^^t events that passed the full re- 
construction chain and all selection criteria used in data. 
The fitted lifetime is consistent with the generated value. 
The simulated events are also used to test that the mea- 
sured proper decay length is not biased with respect to 
the generated one, and that the applied selections have 
the same efficiency for different values of lifetime. 

To test for any bias produced by the fitting procedure, 
500 fast, parameterized Monte Carlo samples are gener- 
ated and analyzed. The average lifetime agrees with the 
generated one, and the assigned uncertainty corresponds 
to the statistical spread of fitted values. 

Another test consists of splitting the data sample into 
two roughly equal parts using various criteria and mea- 
suring the Aj] lifetime in each sample independently. The 
sample is split according to the muon charge, the muon 
direction, the decay length of Kg or the chronological 
date of data taking. All such tests give statistically con- 
sistent values of the AjJ lifetime. 

In conclusion, our measurement of the Aj] lifetime us- 
ing the scmileptonic decay A° — > jiDA'^X results in 
r(AO) = 1.290t^;"|^ (stat)!" ^! (syst) ps. It is consis- 
tent with the current world average lifetime and with 
our measurement in the exclusive decay A^ J/tpA 
[l^ . The D0 results are statistically independent and 
the correlation of systematics between them is very small. 
Their combination results in r(A°) = 1.25llg'|!)gg ps. Our 
new measurements are less consistent with the recent dis- 
crepant measured A^ lifetime than with the current 
world average [J]. 

We thank the staffs at Fermilab and collaborating in- 
stitutions, and acknowledge support from the DOE and 
NSF (USA); CEA and CNRS/IN2P3 (France); FASI, 
Rosatom and RFBR (Russia); CAPES, CNPq, PAPER J, 
FAPESP and FUNDUNESP (Brazil); DAE and DST 
(India); Colciencias (Colombia); CONACyT (Mexico); 
KRF and KOSEF (Korea); CONICET and UBACyT 
(Argentina) ; FOM (The Netherlands) ; Science and Tech- 



nology Facilities Council (United Kingdom); MSMT and 
TABLE 11: Systematic uncertainties in t(A5^) 



Source 


Uncertainty in t(A'[J) 


Detector alignment 


±0.018 ps 


Mass-fitting method 


±0.067 ps 


if-factor determination 


±0.036 ps 


Peaking background 


±0.012 ps 


Resolution scale factor 


±0.036 ps 


Contribution of _B ^ fii/AcX 


+0.000 
-0.027 P° 


Total 


+0.087 
-0.091 P° 



GACR (Czech Republic); CRC Program, CFI, NSERC 
and WestGrid Project (Canada); BMBF and DEC (Ger- 
many); SFI (Ireland); The Swedish Research Council 
(Sweden); CAS and CNSE (China); Alexander von Hum- 
boldt Foundation; and the Marie Curie Program. 



[§: 
b: 

[#: 
[1 
[2: 
[4: 
[s: 
[6: 
[7: 

[8 

[9: 

[10 

[11 

[12 

[13 
[14 

[15 

[16 



Visitor from Augustana College, Sioux Falls, SD, USA. 
Visitor from The University of Liverpool, Liverpool, UK. 
Visitor from ICN-UNAM, Mexico City, Mexico. 
Visitor from Helsinki Institute of Physics, Helsinki, Fin- 
land. 

Visitor from Universitat Ziirich, Ziirich, Switzerland. 

M. Neubert and C.T. Sachrajda, Nucl. Phys. B483, 339 
(1997); M. DiPierro et al, Phys. Lett. B 468, 143 (1999). 

E. Franco et al, Nucl. Phys. B633, 212 (2002). 

F. Gabbiani et al, Phys. Rev. D70, 094031 (2004). 
W.-M. Yao et al. (Particle Data Group), J. Phys. G 33, 
1 (2006). 

CDF Collaboration, A.Abulencia et al, Phys. Rev. Lett. 
98, 122001 (2007). 

DO Collaboration, V. Abazov et al., Nucl. lustrum. Meth- 
ods A 565, 463 (2006). 

V. Abazov et al, Nucl. lustrum. Methods A 552, 372 
(2005). 

D.J. Lange, Nucl. lustrum. Methods A 462, 152 (2001). 
T. Sjostrand et al, Comp. Phys. Commun. 135, 238 
(2001). 

CERN Program Library Long Writeup W5013 (1993); 
documentation available at 

http : // wwwasd . web . cern . ch/ wwwasd/geant / | 

S. Catani, Yu.L. Dokshitzer, M. Olsson, G. Turnock, 

B.R. Webber, Phys. Lett. B269, 432 (1991). 

DELPHI Collaboration, J. Abdallah et al, Eur.Phys.J. 

C32, 185 (2004). 

G. Borisov, Nucl. lustrum. Methods. A 417, 384 (1998). 
DO Collaboration, V. Abazov et al., Phys. Rev. Lett. 94, 
182001 (2005). 

DO Collaboration, V. Abazov et al., Phys. Rev. Lett. 97, 
021802 (2006). 

DO Collaboration, V. Abazov et al., arXiv:hep- 
ex/0704.3909, submitted to Phys. Rev. Lett. 



